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ABSTRACT 



We present new CRIRES spectroscopic observations of the Bry emission line in the nuclear region of the Circinus galaxy, obtained 
with the aim of measuring the black hole (BH) mass with the spectroastrometric technique. The Circinus galaxy is an ideal benchmark 
for the spectroastrometric technique given its proximity and secure BH measurement obtained with the observation of its nuclear 
H 2 maser disk. The kinematical data have been analyzed both with the classical method based on the analysis of the rotation 
curves and with the new method developed by us and based on spectroastrometry. The classical method indicates that the gas disk 
rotates in a gravitational potential resulting from an extended stellar mass distribution and a spatially unresolved dynamical mass 
of (1.7 ± 0.2) x 10 7 Mq, concentrated within r < 7 pc, corresponding to the seeing-limited resolution of the observations. The new 
method is capable of probing the gas rotation at scales which are a factor ~ 3.5 smaller than those probed by the rotation curve analysis, 
highlighting the potential of spectroastrometry. The dynamical mass which is spatially unresolved with the spectroastrometric method 
is a factor ~ 2 smaller, 7.9+] 4 x 10 6 M G indicating that spectroastrometry has been able to spatially resolve the nuclear mass distribution 
down to 2 pc scales. This unresolved mass is still a factor ~ 4.5 larger than the BH mass measurement obtained with the H2O maser 
emission indicating that, even with spectroastrometry, it has not been possible to resolve the sphere of influence of the BH. Based 
on literature data, this spatially unresolved dynamical mass distribution is likely dominated by warm molecular gas and it has been 
tentatively identified with the circum-nuclear torus which prevents a direct view of the central BH in Circinus. This mass distribution, 
■ 1 with a size of ~ 2pc, is similar in shape to that of the star cluster of the Milky Way suggesting that a molecular torus, forming stars at 

, a high rate, might be the earlier evolutionary stage of the nuclear star clusters which are common in late type spirals. 

Key words. Techniques: high angular resolution - Techniques: spectroscopic - Galaxies: active - Galaxies: individual: Circinus - 
Galaxies: kinematics and dynamics - Galaxies:nuclei 



,— ' 1 . Introduction host spher oids like mass, luminosity and stellar velocity disper- 

„ „ , „ , , , , „ sion (e.g. iKor mendv & Richston el Il995l iG ebhardt et al. 2000, 

One of the fundamental yet unresolved problems of modern as- Ferrarese & Merritt 2000 Marconi & Hunt 2003 HLring & mx 

k> . trophysics relates to the formation and evolution of the com- kopjlFerrarese & Ford 2005, Graham 2008a). Moreover, while 

^ ; plex structures that characterize the present-day universe such i{ has , been wMel ted mat Active Galactic Nuclei 

Ctf. as galaxies and clusters of galaxies. Understanding how galaxies (AGN) m ed fe accretion of matter onto a 

formed and how they become the complex systems we observe massive gjj it has recentl been ible to show that BH 

today is therefore a current major theoretical and observational grQWth ifj mosdy due tQ accretion of matter during AGN ac _ 

e 01 ' tivity, and therefore that most ga laxies went through a phase 

There is now strong evidence for the existence of a connec- of strong nudear actiyity (Sokan 19g2 Yu & Xremaine 2002 

tion between supermassive black holes (hereafter BHs), nuclear Marconi et al 2004) It is bdieved that the physical mecha . 

activity and galaxy evolution. Such evidence, that reveals the nism responsible for this coeV olution of BHs and their host 

so-called co-evolution of black holes and their host galaxies, is galaxiefj ifj probably fee dback by the AGN , i.e. the accreting 

provided by the discovery of "relic" BHs in the center of most BH on the host galaxy (SOk & Regs 199g Fabian & Iwasawa 

nearby galaxies, and by the tight scaling relations between BH 1999 Granat0 et al 2004 Di M a tteo et al.ll2005L iMencil 12006. 

mass (M BH ~ 10 6 - 1O IO M ) and the structural parameters of the | Bo wer et alJ^OOrF As the scaling relations between BH mass 

and host galaxy properties represent the clearest sign of co- 
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. „ , , , „ , evolution, it is important to secure these relations by increasing 
Based on observations made with ESO Telescopes at the Paranal 



Observatory under programme ID 383.B-0267(A) 



the number, accuracy and mass range of Mbh measurements. 
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Supermassive BHs are detected and their masses measured 
by studying the kinematics of gas or stars in galaxy nuclei. 
Currently, there are about ~ 50 BH mass measure ments most 
of wh ich are in the ~ 10 7 - 10 9 Mq range (e.g. San i et all 
I2010I) . The majority of these measurements have been made 
with longslit spectroscopy, but the development in recent years 
of Integral Field Unit (hereafter IFU) spectrographs has re- 
sulte d in an improvement in accuracy and reliability (see, 
e.g.. iDavies et alj|2006l iNowak et alj|2007l iNowak et alJi 20 10, 
Krainovic et alJ 120071, iRrainovic et al.| 20091 ICappellar i et alJ 
2009, lNeumaver et al.l2010llRusli et al.l2010l) of the mass deter- 
minations. One crucial issue in BH mass measurements is spa- 
tial resolution: a necessary but not sufficient condition is that the 
spatial resolution must be high enough to make it possible to spa- 
tially resolve the central regions of the galaxy, where the gravi- 
tational effects of the BH dominate those of the host galaxy (the 
so-called 'sphere of influence of the black hole'). Even with the 
advent of Adaptive Optics (AO) assisted observations the best 
spatial resolutions achievable are ~ 0.1", which corresponds to 
~ 10 pc at a distance of 20 Mpc. 

We have developed a new method based on the technique 
of spectroastrometry which allows gas kinematical BH mass 
measurements that partly overcome the spatial resolution limita- 
tions of the "classical" gas (or stellar) kinematical methods. This 
method has been presented and discussed in detail by Gnerucci 
et al. (2010, 2011, hereafter G10, Gil), and provides a simple 
but accurate way to estimate BH masses using either longslit or 
IFU spectra. 

In this paper we apply this spectroastrometric method in an 
attempt to estimate the mass of the nuclear BH in the Circinus 
galaxy and to test whether it is possible to resolve the BH 
sphere of influence. At a distance of 4.2 Mpc ( hence 1" ^ 20 
pc), t his is one of the closest Seyfert 2 galaxies (Freema n et al.l 
1 1977b. Its nuclear activity is indicated by the observed line ratios 
(Olivaetal. 1994), by the ionizatio n cone observed in [O III] 
iMarconi et al.lll994l) and [Si VIII (iPrieto et al.l 120051). by th e 
narrow coronal lines (lOliva et all Il994t Maiolino et al. 2000), 
and by the br oad (FWHM > 3300 km/s) Ha emission in po- 
larized light dOhva et al.lll998l) . The discovery of H2O maser 
emission, whose kinematics are indicative of a circularly rotat- 
ing nuclear disk, has facilitated one of the most impressive and 
accurate mass measurements (1.7+0. 3)xl0 6 Mo ( Greenhil l et al.l 
2003blah of a supermassive BH to date. The relatively low BH 
mass makes the gravitational sphere of influence of the BH very 
small and spatially unresolved even with AO observations at 8m- 
class telescopes. These two facts combined make this galaxy a 
very interesting benchmark to test the potential of the spectroas- 
trometric approach. 

Here we present new long-slit spectra of the nuclear region 
of Circinus obtained with CRIRES at the ESO VLT. This instru- 
ment provides a very high spectral resolution that is essential 
to disentangle the relatively low gas v elocities expected due t o 
the relatively low BH mass estimated by Greenhi ll et al] (12003a). 
We perform two parallel analyses and model the CRIRES spec- 
tra using both the classical rotation curves approach and the 
spectroastrometric method, enabling us to check the quality of 
the modeling and to put tight constraints on Mbh based on a 
comparison of the independent results. In Sect. 2 we briefly re- 
view the key results of G10 and Gl 1 on the application of spec- 
troastrometry to rotating gas disks to enable the detection of nu- 
clear BHs. In Sect. 3 we present the new CRIRES observations 
together with details of the data reduction and subsequent anal- 
ysis. In sect. 4 we describe the CRIRES Circinus gas rotation 
curves and their modeling. In Sect. 5 we present the matching 



spectroastrometric analysis and modeling of the CRIRES spec- 
tra. We compare and discuss the results from the spectroastro- 
metric and rotation curves analyses in Sect. 6, with the aim of 
mutually constraining the two approaches. Finally, in Sect. 7 we 
discuss the results and draw our conclusions. 



2. Spectroastrometric measurements of black holes 
masses 



The spectroastrometric method (see iBailevi Il998l) consists in 
measuring the photocenter of emission lines in different wave- 
length or velocity channels. It has been used by several authors 
to study pre-main sequence binaries and inflows, outflows or the 
disk structure of the gas surrounding pre-main sequence stars 
dTakami et alJl2003HBaines et al.f2 004: Por ter et al.ll2004ll2005L 

IWhelan et aUl2005l). 

In iGnerucci et al. I I20I0I (G10) we illustrated how the tech- 
nique of spectroastrometry can be used to measure black hole 
masses at the center of galaxies. In that paper we focused on 
explaining the basis of the spectroastrometric approach and pre- 
sented an extended and detailed set of simulations showing how 
this method can be used to probe the principal dynamical pa- 
rameters of a nuclear gas disk. A key aspect of G10 was the 
comparison of this technique with the standard method for gas 
kinematical studies based on the gas rotation curve. We demon- 
strated that the two methods are complementary approaches to 
the analysis of spectral data (i.e. the former measures mean posi- 
tions for given spectrum velocity channels while the latter mea- 
sures mean velocities for given slit position channels). The prin- 
cipal limit of the rotation curves method resides in the ability to 
spatially resolve the region where the gravitational potential of 
the BH dominates with respect to the contribution of the stars. 
In G10 we demonstrated that spectroastrometry has the ability 
to provide information on the galaxy gravitational potential at 
scales significantly smaller than the spatial resolution of the ob- 
servations (potentially reaching ~ 1/10). This fundamental fea- 
ture can be illustrated by the following simple example: consider 
two point-like sources located at a distance significantly smaller 
than the spatial resolution of the observations; these sources will 
be seen as spatially unresolved with their relative distance not 
measurable from a conventional image. However, if spectral fea- 
tures, such as absorption or emission lines are present at different 
wavelengths in the spectra of the two sources, the spatial profiles 
of the light distribution extracted at these wavelengths will have 
different centroids, revealing the presence of the two sources. 
The separation between the sources can be estimated from the 
difference in the centroid positions at the wavelengths of the dif- 
ferent spectral features, even if the separation is much smaller 
that the spatial resolution. Thus, the "spectral separation" of the 
two sources makes it possible to overcome the spatial resolution 
limit. 

For clarity, we summarize here the principal features and the 
main steps of the method presented and discussed extensively in 
G10. 

- From the longslit spectrum of a continuum-subtracted emis- 
sion line one constructs the "spectroastrometric curve" by 
measuring the line centroid along the slit for all wavelength 
channels. The "spectroastrometric curve" of the line is given 
by the position centroids as a function of wavelength. 

- From the simulations presented in G10, we showed that 
the information about the BH gravitational field is predom- 
inantly encoded in the "high velocity" (hereafter HV) range 
of the spectroastrometric curve which comprises the points 
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in the red and blue wings of the line. In the case of Keplerian 
rotation, the HV part of the line spectrum originates from gas 
moving at high velocities closer to the BH and this emission 
is spatially unresolved. For this reason the HV part of the line 
spectrum is not strongly influenced by the spatial resolution 
or other instrumental effects like slit losses or by the intrin- 
sic line flux distribution. On the other hand, "low velocity" 
emission (hereafter LV) is usually spatially resolved (i.e. the 
gas moving at lower velocities is located farther away from 
the BH), so the spectroastrometric curve in this regime is 
not useful for the purpose of measuring the BH mass. The 
HV range of the curve is identified by measuring the spa- 
tial extent of the line emission as function of wavelength and 
discarding the central (i.e. LV range) bins where the spatial 
profile becomes broader than the instrumental spatial resolu- 
tion. 

- By measuring the spectroastrometric curves of a given line 
from at least three spectra taken at different slit position an- 
gles, one can obtain a "spectroastrometric map" of the source 
on the plane of the sky by geometrically combining the three 
curves. In the case of integral field spectra, this step is obvi- 
ously not necessary because the spectroastrometric map is 
derived directly from the data cubes. In the case of a ro- 
tating disk with a radially symmetric line flux distribution 
the points of the spectroastrometric map should lie on the 
disk line of nodes. However even for IFU data the effect of 
slit losses or a non-symmetric line flux distribution can per- 
turb the light centroid positions moving them away from the 
disk line of nodes. As shown in Paper I and discussed above, 
these effects become negligible for the HV range of the map 
where the emission is spatially unresolved. In contrast, the 
LV points of the map tend to lie away from the line of nodes 
in a typical "loop" shape. The final spectroastrometric map 
is then obtained by selecting only the HV points. 

- One can then estimate the disk line of nodes by a linear fit 
to the HV range of the spectroastrometric map, project those 
points on the estimated disk line of nodes and obtain the disk 
rotation curve. Finally, one can apply a simple model fitting 
procedure to obtain the parameters determining the gas rota- 
tion curve, and in particular the BH mass. 

In a subsequent work (Gnerucci et al. 20 1 1 , G 1 1 ), we applied 
the method to the radio galaxy Centaurus A, using seeing lim- 
ited longslit spectra and both seeing limited and AO assisted IFU 
spectra. We selected this galaxy because it has been extensively 
studied using the gas kinematical method; the observed kinemat- 
ics show that the gas is circularly rotating around the BH and the 
mass and the disk parameters are well constrained. Moreover 
both longslit and IFU data with different spatial resolutions are 
available, allowing a direct comparison of the spectroastrome- 
try results to the results obtained from data sets with different 
characteristics. We demonstrated with real data that the method 
provides the same mass estimates as the classical gas kinemati- 
cal method, but allows smaller spatial scales to be probed . 

We briefly summarize the key results of Gl 1 here. 

- We demonstrated with real data the capability of spectroas- 
trometry to overcome the spatial resolution limit. The mini- 
mum distance from the BH at which we can measure the gas 
rotation curve is ~ 20mas, which corresponds to ~ 1/15 of 
the spatial resolution for the seeing limited data and ~ 1/6 
for the AO assisted data. 

- The method based on spectroastrometry is in excellent agree- 
ment with the classical method based on the rotation curves 
(the BH mass estimates are consistent within +0.2 dex and 
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Fig. 1. Continuum subtracted detector 3 CRIRES spectrum (slit 
PA1) extracted at the position of the continuum peak. Solid red 
line: fit to Bry using a Gauss-Hermite function. 

within less than +0.05 dex when comparing the results ob- 
tained from the same dataset). This demonstrates the robust- 
ness of the method: using the same dataset we can indepen- 
dently apply the classical and the spectroastrometric meth- 
ods and obtain consistent results. 

- The application of the spectroastrometric method to differ- 
ent types of data (IFU and longslit) give consistent results 
(within only ~ 0. 1 dex), demonstrating the versatility of our 
method. 

- A typical feature of spectroastrometry is its insensitivity to 
disk inclination. In our modeling, in fact, mass and disk in- 
clination are coupled. Therefore our method can return only 
a value for Mbh sin 2 i; the inclination must be assumed or de- 
termined by another method to derive the value of the mass. 

- The spectroastrometric method has been applied indepen- 
dently of the classical method. However as discussed in G10, 
it is clear that the spectroastrometric and "classical" rotation 
curves are complementary and orthogonal descriptions of the 
position velocity diagram. Therefore, a development of this 
method will be its application in combination with the clas- 
sical method based on rotation curves. This provides con- 
straints on the disk inclination and thus removes the mass- 
inclination degeneracy. 

The utility of the method is not limited to gas disks rotating 
around supermassive black holes but also to galaxy disks and 
indeed we have applied the method to improve virial mass es- 
timates of high-z galaxies obse rved with IFU spectra (for more 
details see lGnerucci et alfcoi lb V 

3. Observations and data reduction 

Circinus was observed between April and July 2009 (pro- 
gram 383.B-0267A) with CRIR ES (Cryogenic hig h-resolution 
Infrared Echelle Spectrograph ; iKaeufl et al.l 120041) at the ESO 
VLT with the aid of MACAO (Multi Applications Curvature 
Adaptive optics) Adaptive Optic system. The observations con- 
sist of three nuclear spectra characterized by a different slit po- 
sition angle: the slits with position angles of 30°, 90° and 150° 
with respect to the north direction are denoted by PA1, PA2 and 
PA3, respectively. The orientation of PA1 was chosen in order to 
align the slit with the major axis of the galaxy, and PA2 and PA3 
were incrementally rotated by 60° in order to provide a uniform 
coverage in position angles with the available slits. 

The Adaptive Optics module MACAO was used with the 
unresolved Circinus nucleus as a guide star. The spectra were 
obtained with a 0.4" wide slit (slit length ~ 31") and a re- 
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Fig. 2. Rotation curves for the CRIRES Bry spectra. Left panel: slit "PA1". Central panel: slit "PA2". Right panel: slit "PA3". The 
solid red lines represent the curves expected from the model. 



solving power of A/AA = 50000, corresponding to a spectral 
resolution of ~ 0.4A (~ 5.5km/s) at the observed Bry wave- 
length (2.169fim). The resulting wavelength range is ~ 2.137 - 
2.188 fim, covered with a mosaic of four InSb arrays providing 
an effective 4096 x 512 focal plane detector array in the focal 
plane. The choice of the very high-spectral resolution is criti- 
cal for the measurement of the small velocities and velocity dis- 
persions which are exp ected, given the relatively low BH mass 
dGreenhill et alj2003al) . The detector pixel scale is 0.086"/pixel 
along the slit axis and 0.102A/pixel along the dispersion direc- 
tion. For each position angle we obtained 6 spectra with an expo- 
sure of 180s, nodding the target along the slit for sky subtraction, 
resulting in a total on-source exposure of 1080s. 

The data were reduced using the CRIRES reduction pipeline 
(version 1.8) for the subtraction of dark frames, flat-field correc- 
tion, sky subtraction, wavelength calibration and for the combi- 
nation of the various observing blocks. IRAF tasks3 were then 
used to rectify the spectra and to correct for telluric absorption 
using a suitable standard star. During data reduction, we con- 
centrated only on the third detector that targeted the Bry line 
(rest frame wavelength 2.1661yura) observed at ~ 2.169^/m. The 
spatial resolution of the data was estimated in the final reduced 
spectra from the extension of the nuclear continuum along the 
slit (the galaxy nucleus with an estimated size of 1 .9 + 0.6 pc, 
corresponding to 071 + 0703 at 3.5 Mpc, is i n fact unresolved a t 
the spatial resolution of these observations: iPrieto et al.l 12004). 
We obtain an average value of ~ 0.7" FWHM for the spectra at 
the three PAs indicating that the AO correction did not perform 
well, probably due to the non point-like nature of the Circinus 
nucleus in the R band. The spectra were continuum subtracted 
since, as discussed in G10 and Gil, the spectroastrometric anal- 
ysis requires a continuum subtracted spectrum. Finally, as the 
dispersion was higher than needed (~ 1 .4 km/s/pixel at Bry), we 
rebinned the spectra by 5 pixels along the dispersion direction to 
improve the signal-to-noise ratio. 

The final continuum subtracted Bry spectra have a typical 
signal to noise ratio (hereafter S/N) ~ 35. This value is obtained 
from the nuclear spectrum extracted from an aperture along the 
slit equal to the spatial resolution; the Bry profile is fitted with a 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



Gauss-Hermite model (see below) and the quoted S/N is defined 
by the ratio between the peak model value and the rms of the fit 
residuals. 

With our observational setup (10 km s _1 spectral resolution, 
077 spatial resolution) the minimum detectable BH mass is de- 
termined by the spatial resolution and the thermal broadening of 
the ionized gas (i.e. <x ~ 15-22 km for T e = (1 — 2) X 10 4 
K), and corresponds to 

MBH = ^- L9Xl ° 5M0 (2Ok^) 2 (o^) (1) 

if r = 072 is the spatial scale probed with spectroastrometry. At 
this stage, this is a purely indicative value. 

4. Analysis of the rotation curves 

We measure the gas rotation curves in the three spectra (PA1, 
PA2 and PA3) by fitting the Bry line r ow by row along the slit 
direct ion with a Gauss-Hermite profile (Cappell ari & Emselleml 
2004) thus obtaining mean velocity, velocity dispersion and line 
flux at each position along the slit. In Fig. Q] we show a fit exam- 
ple corresponding to the "PA1" continuum subtracted spectrum, 
rebinned by 5 pixels along the dispersion direction and extracted 
from one pixel along the slit at the position of the continuum 
peak. The rotation curves obtained at the three slit position an- 
gles are presented in Fig. [2] and clearly show the typical S-like 
shapes expected from rotating disks. Fig.[3]shows the Bry sur- 
face brightness profiles along the slits (black lines) compared 
with the corresponding continuum curves (red lines). Since con- 
tinuum emission is dominated by the unresolved ho t dust emis- 
sion at the spatial resolution of our observations dPrieto et al.l 
2004), figures [3] and [2] clearly shows that we spatially resolve 
the emission regions and the kinematics of Bry . 

The rotation curves are first analyzed according to the "clas- 
sical" gas kinematical method. The modeling of the rotation 
curves is performed fo llowing the procedure first described by 
iMacchetto et al.l (1 19971) and subsequently refined by several au- 
thors; here we used the m odeling code described in detail by 
Marc oni et al.l (2003. 2006). Very briefly, the code computes the 
rotation curves of the gas assuming that the gas is rotating in 
circular orbits within a thin disk in the principal plane of the 
galaxy potential. We neglect any hydrodynamical effect like gas 
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Fig. 3. Surface brightness profiles of the Bry line along the slits 
(black lines). The red lines represent the surface brightness of 
the underlying continuum. 



pressure. The gravitational potential is made up of two com- 
ponents: the stellar potential, characterized by its mass-to-light 
ratio M/L, and a spatially unresolved mass concentration char- 
acterized by its total mass Mjjdm (UDM means "Unresolved 
Dynamical Mass"). This spatially unresolved dynamical mass 
concentration is dark in the sense that it is not accounted for 
by the continuum emission used to trace the stellar mass distri- 
bution. It could represent the putative supermassive black hole 
but also any spatially unresolved stellar or gaseous component 
or a combination of them. In computing the rotation curves we 
take into account the finite spatial resolution of the observations, 
the intrinsic surface brightness distribution of the emission lines 
(hereafter ISBD) and we integrate over the slit and pixel area. 
The free parameters characterizing the best fitting model are 
found by standard^ 2 minimization. 

We first determined the stellar gravitational potential from 
a K-band surface brightness profile of the Circinus Galaxy 
which was obtained from a NICMOS/CAMERA2 F222M mo- 
saic kindly provided by Roberto Maiolino (for detai ls on ob- 
servations and data reduction see Mai olino et al.ll2000l) . The in- 
version procedure to derive the stellar density distribution from 
the surface brightness is not unique if the gravitational poten- 
tial does not have spherical symmetry. Assuming that the grav- 
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Fig. 4. Upper panel. Fit of the NICMOS F222m light profile. 
The black continuous line represent the best-fit model. The dot- 
dashed line represents the profile of the compact component 
(0715 FWHM size) and the dashed line the extended compo- 
nent. The compact component has been detected a nd identified 
with h ot dust emission in NACO observations by iPrieto et all 
(2004). Note that the spatial extension of the compact compo- 
nent is similar to that of the NICMOS PSF and therefore the first 
Airy ring is still visible. Lower panel: model residuals. 



itational potential is an oblate spheroid, the inversion requires 
knowledge of the potential axial ratio q, and the inclination of 
its principal plane with respect to the line of sight, i. These 
two quantities are related by the observed isophote ellipticity, 
and for the Circinus galaxy we solved this ambiguity by im- 
posing an axial ratio of ~ 0.1 that corresponds to an almost 
disk-like spheroid; Circinus is in fact an SA(s)b spiral galaxy 
and, moreover, the stellar contribution to the gravitational poten- 
tial is maximized by imposing a disk-l ike geometry. Following 
Ivan der Marel & van den Bosch! d!998l) . we then converted the 
oblate spheroid density distribution model to an observed sur- 
face brightness distribution in the plane of the sky by integrat- 
ing along the line of sight, convolving with the Point Spread 
Function (PSF) of the telescope+instrument system and averag- 
ing over the detector pixel size. Then, the derived model light 
profile can be directly compared with the observed one. A de- 
tailed description of the relev ant formulae a n d the inversion and 
fit procedure is presented in iMarconi et al ] (I2003h . The results 
of the analysis of the NICMOS F222M profile are presented in 
Fig. |4] In the upper panel the black solid line represents the 
model surface brightness profile convolved for the instrumental 
effects. This model profile is the combination of an extended and 
a compact components, represented by dashed and dot-dashed 
lines, respectively. The extended component represents stellar 
emission and is parameterized with the oblate spheroid described 
above. The compact component is characterized by spherical 
symmetry and by a sharp boundary, resulting in a FWHM of 
0/15. This component is consistent with th e compact K band 
source of size ~ 0/2 detected with NACO bv lPrieto etall (f2004) 
and identified with hot dust emission from the AGN torus on the 
basis of the steep near-IR colors. This identificati on has been fur- 
ther confirmed by Muel ler Sanchez et al .l J2006h who were able 
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Fig. 5.x 1 vs. i curve for the determination of best value and con- 
fidence interval of the disk inclination. The dotted lines repre- 
sents the 1, 2 and 3cr confidence levels. 

to deconvolve the stellar and dust continua, confirming that the 
compact K band emission is dominated by hot dust. Therefore 
we have not considered it for the stellar mass profile used in the 
analysis of the rotation curves. 

The best fitting rotating disk model is represented in Fig. [2] 
by the red solid line and, overall, it nicely reproduces the ob- 
served rotation curves. One possible cause of concern is that the 
rotation curve at PA3 shows a ~ 50km/ s velocity shift with re- 
spect to the rotating disk model and, indeed, the red line in the 
right panel of Fig. [2]has been offset by the same amount to match 
the observed points. PA3, with a position angle of 150° is almost 
perpendicular to the line of nodes, which approximately aligns 
with the PA1 slit (P.A. 30°). This shift might then be accounted 
for by outflowing motions since the PA3 slit is lying within the 
ionization cone illuminated by the central AGN. Such outflows 
are commonly observed in Seyfert galaxi es and are likely driven 
by AG N radiation pressure or winds (e.g. lMueller-Sanchez et alj 
1201 ll and references therein). However, since the rotation curves 
from the PA1 and PA2 slits are both well reproduced by a rotat- 
ing disk model and since, after allowing for the velocity shift, 
the PA3 slit rotation curve is also well reproduced by the same 
model, we believe that the outflowing motions detected in PA3 
do not affect the fin al mass estimate, simil arly to what has been 
found in Cygnus A dTadhunter et al.ll2003l) . 

Since the disk inclination is coupled with the unresolved 
dynamical mass (i. e. Mudm and stellar M/L ratio, see, e.g., 
iMarconi et al.ll2006h . we choose to keep it as a fixed parame- 
ter in the modeling to avoid convergence problems with the x 1 
minimization. Its best value and confidence interval is then iden- 
tified by using a grid of fixed i values and performing the fit for 
each of them. In this way we construct the^ 2 vs. i curve shown 
in Fig. [5] which allows us to identify both the be st i value and 
the confidence intervals (for details see lAvnill 19761) . The best in- 
clination value and confidence interval (at lcr confidence level) 
is (61 + 1)°. Thi s value is in good agre ement with the value of 
65° recovered bv lFreeman et alj d 19771) from large- scale optical 
photom etry and with the value of 55° obtained by iHicks et alj 
(2009) for the inner H2 gas disk (radius < 9pc). 

The best-fit model parameters for the Bry rotation curves are 
reported in Table. [1] while in Fig. [2] we show the Bry rotation 
curves and the relative best-fit model. We estimate the uncertain- 
ties on the best-fit model parameters using the bootstrap method 
dEfron & Tibshirani 1994). The dataset consists of 181 elements 
in total from the rotation curves of the three slits, each character- 
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Fig. 6. Confidence regions in the Mudm-M/L plane for the ro- 
tation curve model fits. The iso-contours indicate the 95%, 90% 
and 68% (the inner one) confidence levels. The dotted line marks 
the Mudm-M/L values corresponding to the minimum;^ 2 fit. The 
open circles indicate the grid of values used in the Mudm-M/L 
plane. 



ized by spatial position, flux, velocity and velocity dispersion. 
We randomly extracted from this dataset a subsample with the 
same number of elements. Due to the random extraction, each 
subsample will have some elements replicated a few times and 
some elements that are entirely missing. We then performed the 
fit on these randomly extracted 100 subsamples and estimated 
errors on parameters by taking the standard deviation of the best 
fit values, usually normally distributed. 

Since Mudm and M/L are correlated, we have also evalu- 
ated the confidence intervals for the coupled distribution of these 
two model parameters. To accomplish this we constructed a grid 
of values in the Mudm - M/L plane and performed the fit for 
each of these points by varying all the other parameters. In this 
way we created a x 1 surface which allows us to identify both 
the best values and the confiden ce regions in the Mudm - M/L 
plane(for details see lAvnill 19761) . In Fig.|6]we show the resulting 
contours of the confidence regions and each contour is labeled 
with the corresponding confidence level. As shown in the figure, 
the confidence regions are consistent with the statistical uncer- 
tainties obtained with the bootstrap method for these parameters 
(cf. Table [TJ and the Mudm and M/L parameters show the ex- 
pected anti-correlation. 

5. Spectroastrometric analysis 

5. 1 . Analysis of the spectra 

As discussed in G10 and Gil to reduce the effects of noise in 
measuring the spectroastrometric curves, we constructed a "syn- 
thetic" position velocity diagram (hereafter PVD) of the Bry 
line for each longslit spectrum. We fitted the Bry line with a 
Gauss-Hermite function at all positions along each the PA1, PA2 
and PA3 slits, obtaining the relative model parameters (flux, 
mean velocity, velocity dispersion and Hermite parameters /13 
and /14). We then use the fitted profiles of Bry to reconstruct 
noise-free synthetic PVD's for each slit. 
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The PVD for the original Bry longslit spectrum at PA1 is 
compared with its "synthetic" version in Fig. [7] It can be seen 
that the "synthetic" PVD is noise-free because each row rep- 
resents the fitted parametric profile. However, one has to take 
into account the errors on the fitted parameters in order to es- 
timate the errors on the synthetic pixel counts. Therefore, for 
each spectral fit, we simulated 1000 synthetic spectra from 1000 
realizations of the set of the five profile parameters distributed 
following a pentavariate distribution with the fit correlation ma- 
trix. The flux and error of each pixel is then estimated from the 
mean and standard deviation of the 1000 realizations. We note 
that for each flux profile along the slit, the errors on fluxes are 
uncorrelated because they originate from independent fits to the 
line profiles. 
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Fig. 7. Position velocity diagram of the continuum subtracted 
Bry spectrum at PA1. Upper panel: observed. Bottom panel: 
"synthetic" reconstruction. The horizontal dotted line overplot- 
ted on each panel represents the continuum peak position. The 
isophotes denote the same values in both panels. 

From the "synthetic" Bry spectra for the three slits (PAl, 
PAl and PA3), obtained following the method outlined above, 
one can derive the spectroastrometric curves which are shown in 
Kg.! 

From Fig. [8] we can see that the position of the light cen- 
troids in the high velocity (HV) range tend to approach the po- 
sition (i.e. the position of the continuum peak). As discussed in 
G10 and Gil, this is the expected behavior of spectroastromet- 
ric curves in the presence of an unresolved central mass (i.e. the 
BH). This behavior is less evident on the red side because of the 
small number of high velocity points. Errors on photocenter po- 
sitions range from ~ 0.02" to ~ 0.1", that is from ~ 1/35 to 
~ 1/5 of the spatial resolution of the data; this is the accuracy 
with which we can measure centroid positions. 



5.2. The spectroastrometric map of the source 

We obtained three spectroastrometric curves from Bry , one for 
each PA of the slit. Each spectroastrometric curve provides the 
photocenter position along one slit, i.e. the position of the pho- 
tocenter projected along the axis defined by the slit direction. 
Thus, combining the spectroastrometric curves obtained at dif- 
ferent PA, we can obtain the map of photocenter positions on 
the plane of the sky for each velocity bin. In principle, the spec- 
troastrometric curves from two non-parallel slits should suffice 
but we can use the redundant information from the three slits to 
recover the 2D sky map as described in Section 4.1 of G10. We 
have chosen a reference frame in the plane of the sky centered 
on the center of the PAl slit (which corresponds to the position 
of the continuum peak along the slit) with the X axis along the 
North direction. For a given velocity bin we then determined the 
position of the light centroid on the sky plane, producing the 2D 
spectroastrometric map shown in Fig. [9] 




ARA (arcsec) 

Fig. 9. 2D spectroastrometric map for the Bry line. The red point 
marks the inferred position of the point-like mass. The red solid 
line represents the line of nodes of the disk obtained from a lin- 
ear fit to the HV points, as described in the text. The dotted red 
lines represent the 1 <x uncertainties on the line of nodes position 
angle. The origin of the map is the position on the sky plane of 
the PAl slit center (which corresponds to the continuum peak 
position). 

As discussed in G10 and Gil the coordinates on the plane 
of the sky of the center of the PA2 and PA3 slits must be con- 
sidered as free parameters. These unknowns are estimated si- 
multaneously with the position of the photocenter following a 
X 2 minimization procedure. The final spectroastrometric map on 
the plane of the sky is that given by the best fitting set of slit 
centers. The error bars on the points represent the uncertainties 
resulting from the fit. The black points correspond to the "high 
velocity" points (i.e. v < 295 km/s and v > 475 km/i) which 
were actually used to determine the location of the slit centers. 
The HV parts of the line profile represent the tails of the line- 
of-sight velocity distribution and could not be selected by iden- 
tifying spatially unresolved emission along the slit, i.e. the cri- 
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Fig. 8. Spectroastrometric curves of the Bry line at the three slit position angles. Left panel: PA1. Central panel: PA2. Right panel: 
PA3. The dashed vertical lines on each panel represent the limits of the "high velocity" range. Wavelengths have been converted 
into velocity using as reference (zero velocity) the rest frame Bry wavelength (2.1661/im). 



terion adopted in Gl 1, since emission along the slit direction is 
spatially resolved. We therefore extracted a spectrum of the Bry 
line by averaging the CRIRES spectrum over an aperture equal 
to the spatial resolution and centered on the continuum peak and 
we selected the velocity bins characterized by a flux lower than 
0.3 times the line peak flux. This criterion to identify the HV 
range is different from the one used in Gl 1 for the Centaurus A 
spectra. However, it produces the same HV range for Centaurus 
A and is therefore equivalent. The spatial extension of the emis- 
sion in the HV region most probably indicates that the high ve- 
locity spectroastrometric measurement probes the gravitational 
potential of an extended mass distribution rather than that of the 
BH. This issue will be discussed below, where we will show that, 
indeed, the spectroastrometric analysis indicates the presence of 
an extended mass distribution. The emission of the gas moving 
at higher velocities does not originate from the inner region of 
the disk rotating around the BH but feels the gravitational poten- 
tial of an extended mass component. 

The two-dimensional spectroastrometric map just described 
(and shown in Fig. [9} can now be used to estimate the geomet- 
rical parameters of the nuclear gas disk. If the gas kinematics 
are dominated by rotation around a point-like mass {Mudm), the 
position of the light centroid at the high velocities should lie on 
a straight line (which identifies the direction of the disk line of 
nodes) and should approach, at increasing velocities, the posi- 
tion of the rotation center. This consideration allows us to make 
a first estimate of the position on the plane of the sky of the ro- 
tation center as the average position of the HV points in the 2d 
spectroastrometric map2 as well as the position angle of the line 
of nodes (9lon), obtained by fitting a straight line to the high 
velocity points (see Fig. A more accurate estimate of these 
parameters, derived from the model fitting, will be described in 
the next section (see TableQ]). The uncertainty on center position 
is ~ 0.04", i.e. ~ 1/18 of the spatial resolution of the data. 



2 In practice, we calculate this position by first taking the averages 
of the coordinates of the points in the blue and red HV ranges and then 
taking the average coordinates of the "blue" and "red" positions. 



5.3. Estimate of the unresolved mass from the 
spectroastrometric map 

We recover the value of the point-like mass from the spectroas- 
trometric map, following the method presented and discussed in 
detail in G10 and Gil. 

Briefly, under the assumption that the gas lies in a thin disk 
configuration inclined by i with respect to the plane of the sky 
(i = face-on) and the disk line of nodes has a position angle 
Own, the line of sight component of the circular velocity of a 
gas particle with distance r from the rotation center (hereafter 
V c h for "channel velocity") is given by: 



Vch 



G[M UDM + MIL ■ L(r)] . 
+ -i / sin(0 + V sys 



(2) 



where L{r) is the radial luminosity density distribution in the 
galactic nucleus, M/L is the mass to light ratio of the stars and 
we also added the systemic velocity of the galaxy V sys . 

First we recover the disk line of nodes by fitting a straight 
line on the 2d map. Then we project the position of the 2d map 
points (x c h,y c h) on the line of nodes, calculating their coordinate 
with respect to this reference axis (S c h) and then their distance r 
from the rotation center used in Eq.[2](ie. r = k\S c ;, - Sol where 
S o is the coordinate along the line of nodes of the disk and k is 
a scale factor to transform angular to linear distance^ see G10 
for details). 

The unknown parameters of this model are found by mini- 
mizing the quantity 



X 



= 1 

ch 



v ch - v ch 



model 



A(S ch ;par) 



(3) 



where A(S c h',par) is the uncertainty of the numerator. As 
previously discussed, we restrict the fit (i.e. the sum over the ve- 
locity channels) to the HV range of the map. The channel veloc- 
ity V c han has no associated uncertainty since it is not a measured 
quantity but the central value of the velocity bin. Finally we add 
a constant error (A SVJ ) in quadrature to the quantity A{S c i,; par) 



3 For consistency with previous works we assume a distance to 
Circinus of 4.2 Mpc. At this distance 1" corresponds to ~ 20Apc. 
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in order to obtain a reduced^ 2 close to 1 (see sect. 5 of G10 for 
a detailed explanation of this choice). Finally, using the best fit 
values of the model parameters we can compute the (x, y) posi- 
tion of the rotation center in the sky plane. 

We verified that the value of Mudm is affected little by the 
actual values of the position angle of the line of nodes and of So, 
the position of the rotation center along the line of nodes. Indeed, 
repeating the fit with 9lon and So values randomly varied within 
the uncertainties, Mudm changes by ~ 0.06 dex, consistently 
with the typical lcr statistical uncertainties of the fit (see Table 

As discussed previously and in G10 and Gil, disk incli- 
nation i and spatially unresolved mass concentration are cou- 
pled and this fitting method can only measure Mjjdm sin 2 i. 
Therefore, we assumed a value for the inclination to obtain a 
mass value. In the following, we will consider a disk inclination 
i = 61°, as obtained in the rotation curve analysis. Summarizing, 
the free parameters of our fit are: 

Mudm spatially unresolved mass; 
M/L mass to light ratio of the nuclear stars; 
So position of the rotation center along the line of 
nodes ; 

Vsys systemic velocity of the galaxy; 

The fit results are tabulated in Table [TJ and presented graph- 
ically in Fig. [10] where we plot the r = \S c h - Sol vs V — 
\V c h - Vsy S \ rotation curve and the line of nodes projected ro- 
tation curve (e.g. S c i, vs V c fd- The solid red lines represent the 
curves expected from the model (r vs \V c i" wdel - V sys \ and S c h vs 

V ch " wdel respectively). 

The first important result from this analysis is that the mini- 
mum distance from the rotation center at which there is a veloc- 
ity estimate is ~ 0.1" corresponding to ~ 2.3 pc (see right panel 
of Fig. [Tol l, while with the standard rotation curve method the 
minimum distance from the center which can be observed is of 
the order of half the spatial resolution (~ 0.35"). This clearly 
shows how spectroastrometry can overcome the spatial resolu- 
tion limit. 

It is not possible to obtain measurements at spatial scales 
smaller than ~ 1 /3 of those probed by the rotation curve analysis 
because of the limited signal-to-noise ratio of the spectra; this 
prevents us from measuring the photocenters of the line at higher 
velocities, further extending the curves in Fig. [8] to the red and 
to the blue sides. 

6. Comparison of the rotation curves and 
spectroastrometric analyses 

Here we compare the results of the two independent analyses of 
the CRIRES spectra (i.e. the fits labeled as "Fit-R 0" and "Fit- 
S 0" in Table [TJ and discuss the combined constraints from the 
two methods. Eventually, we will re-analyze the data taking into 
account these combined constraints. 

As discussed in G10 and Gil, the modeling of the spectroas- 
trometric data is not able to constrain the mass to light ratio M/L 
of the nuclear star distribution because it models the gas rotation 
curves at small distances from the rotation center where the con- 
tribution of the stellar mass to the gravitational potential is neg- 
ligible. In fact performing a first model fit with M/L free to vary 
("Fit-S 0" in Table [TJ we obtain a value of log w (M/L) ^ -11. 
Therefore we used the results of the rotation curve modeling to 
constrain this parameter. In the "Fit-S 1" presented in Fig. [10] 
and TablefJJwe fixed M/L to the value obtained from the rotation 



curves analysis. Another important check concerns the geometri- 
cal parameters of the modeled gas disk (i.e. the position angle of 
the line of nodes Own an d the disk center, or equally, the center 
position on the plane of the sky (x c ,y c )). The spectroastrometric 
model fitting gives a value of 217° for Olon- This is consistent 
with the rotation curve modeling, which gives a best-fit value 
of 215°. On the other hand, the rotation curve modeling gives a 
rotation center that is shifted south by ~ 0.05" with respect to 
that returned by the spectroastrometric modeling. However, the 
spectroastrometric map probes a spatial region of radius ~ 0.2", 
whereas the rotation curves probe distances up to ~ 3" from the 
disk center; therefore, we expect the former to give a more robust 
indication of the center position than the latter. 

A further constraint that the spectroastrometry can place on 
the rotation curves analysis is the position of the centers of the 
three slits on the plane of the sky. As previously described, when 
we combine the three spectroastrometric curves to obtain the 
spectroastrometric map, we obtain an estimate of the positions 
of the slit center with an accuracy better than that of the tele- 
scope pointing (see G10 and Gl 1 for details). 

Therefore, we repeated the rotation curves modeling impos- 
ing geometrical constraints obtained from the spectroastrometric 
analysis, namely the slit center positions and the values of Own 
and (x c ,y c ) ("Fit-R 1" in TablefTJ. In fact, this produces a higher 
value of the reduced^ 2 than the original fit ("Fit-R 0" in Table 
[TJ|. However, the best-fit Mudm and M/L values are the same, 
within the uncertainties. The differences relative to the "Fit-R 

0" case ( 0.03 dex and of +0.03 dex, for M UD m and M/L, 

respectively; table [JJl can be considered an indication of the in- 
fluence of the geometrical disk parameters on the derived Mudm 
and M/L values. 

The slightly worse fit quality of "Fit-R 1" may also be an 
indication that that the inner disk probed by the spectroastro- 
metric analysis has a different geometry (in this particular case 
the disk geometrical center position) than the outer disk, which 
influences the rotation curves. The rotation curves, in fact, probe 
the average disk geometric parameters over a spatial region that 
extends over the central ~ 100 pc whereas spectroastrometry 
probes the inner ~ 10 pc. 

Regarding the principal parameter of the modeling, Mudm, 
we note that the rotation curve fits give a mass estimate larger by 
» 0.35 dex with respect to the values returned by the spectroas- 
trometric modeling. This difference is larger than the statistical 
or systematic uncertainties, but the good agreement between the 
two modeling approaches with respect to the geometrical disk 
parameters gives confidence that the results of the fitting are ro- 
bust. We therefore attribute the difference in Mudm to the dif- 
ferent spatial scales probed by the two methods. For the rotation 
curves approach, as previously noted, the spatial resolution of 
the observations sets the minimum radius at which the gas rota- 
tion curve can be probed. Therefore, with this approach, we can 
only measure the total dynamical mass enclosed within a radius 
equal to half the spatial resolution, in this case, corresponding 
to ~ 0.35" or 7 pc. On the other hand, the spectroastromet- 
ric approach is able to probe the gas rotation curve at smaller 
radii, allowing us to measure the total dynamical mass enclosed 
within ~ 0.1" (~ 2 pc, see Fig.QJJi. We conclude that the differ- 
ence in the mass measured by the two approaches (7.6 x 10 6 M o 
from spectroastrometry and 1.7 x 1O 7 M from rotation curves) 
is most likely due to an extended mass distribution within the 
inner ~ 0.35" (7 pc), which is large compared to the BH mass 
(~ 1.7 x 1O 6 M ). 
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V ch (km/s) r (arcsec) 

Fig. 10. Results of the fit to the spectroastrometric Bry data. The left panels show the line of nodes projected rotation curve S c t, vs. 
V c h- The right panels show the r = \S c /, - Sol vs. V = |V C /, - V sys \ rotation curve. The solid red lines represent the curves expected 
from the model. The dotted line on right panel represents the rotation curve due to the star component of the gravitational potential. 

Table 1. Fit Results. 



Parameter 




Best fit value±error 






Bry Rotation Curves Modeling (1) 


Bry Spectroastrometric Modeling 




Fit-R 


Fit-R 1 


Fit-S 


Fit-S 1 


blon n 


215.0 ±0.5 


217 < 2 > 


217[±4] (3) 


217[±4] (3) 


X C ["] 


-0.013 ±0.02 


0.04 (2) 


0.04 ± 0.04 


0.04 ± 0.03 


y c ["] 


0.067 ± 0.02 


0.08 (2) 


0.08 ± 0.03 <4) 


0.08 ± 0.02 (4) 


log l0 (M UDM /M o ) 


7.22 ± 0.05 


7.19 


6.93 ± 0.07 


6.88 ± 0.07 


log l0 (M/L) 


-0.74 ±0.01 


-0.737 


-11.2 ±0.0 (5) 


-0.74 < 2 > 


i □ 


61[±1] (3) 


61[±1] (3) 


61[±1] (3) 


61[±1] (3) 


V sys [km/s] 


389.3 ± 0.8 


393.6 


391 ± 16 


390 ±9 


XL (X 2 /D.O.F.) 


1.80 (309/172) 


3.69 (645.4/175) 


0.25 (0.76/3) 


0.3 (1.4/4) 



Notes. (1) The best-fit parameter confidence intervals are computed only for the interesting parameters as explained in the text (see also I Avnil 
119761) . (2) Parameter held fixed. (3) Parameter held fixed. The adopted value (with corresponding errors) was estimated independently from the fit 
presented in the table. (4) Derived from x c using the adopted value of 9 L q N . (5) Parameter not constrained from the fit. 



7. The nuclear mass distribution 

The most important result from our attempts to measure the mass 
of the BH in Circinus is that the Mudm value estimated from ro- 
tation curve modeling is ~ 0.35 dex higher than the value es- 
timated from spec troastrometric modelin g. Moreover, the BH 
mass measur ed bv lGreenhilletalJl2003al from the H2O maser 
spectroscopy (1.7x 1O 6 M ) is even lower (~ 0.7-1 dex) than our 
estimates. The analyses based on rotation curves and spectroas- 
trometry can robustly probe gas rotation and thus the mass dis- 
tribution, at different distances from the dynamical center. The 
lower limits on these distances ("limiting inner radii") are set by 
the instrumental setup. In the case of rotation curves the limiting 
inner radius is half the instrumental spatial resolution and corre- 
sponds to ~ 0.35" or 7pc for our Circinus data. In the case of 



the spectroastrometric analysis the limiting inner radius is a frac- 
tion of the spatial resolution and corresponds to ~ 0.1" or 2 pc 
(right panel of Fig.lTOt. Therefore the different measurements of 
the unresolved dynamical mass indicate the existence of an ex- 
tended mass distribution, revealed by the different spatial scales 
probed. This extended mass distribution could not be identified 
in the K-band continuum light profile, also for the presence of 
the st rong hot dust emission within the inner 072 dPrieto et alJ 
2004). It was only revealed as unresolved mass in both our dy- 
namical analyses. 

In Fig.QT]we summarize existing dynamical mass measure- 
ments performed in the nuclear region of the Circinus galaxy. 
The enclosed mass is plotted as a function of distance from the 
dynamical center, which is assumed to coincide with the BH lo- 
cation. The blue point represents the BH mass measurement by 
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Fig. 11. Dynamical mass measurements in the nuclear region 
of the Circinus galaxy. Measurements are marked by filled cir- 
cles of different colors and are identified in the legend. The 
red filled circle represents the CRIRES spectroastrometry mass 
measurement obtained assuming a disk with i = 61° with re- 
spect to the line of sight, same as large scale measurements; 
the open red circle represents the spectroastrometry mass mea- 
surement obtained assuming an edge-on disk, same as the maser 
disk. The region shaded by gray dotted lines is characterized by 
r < 20 mas and is the region probed by the H2O maser rota- 
tion curve, where the dynamical mass is mostly unresolved and 
belongs to the BH. The enclosed mass within this region must 
be dominated by the BH, since no evidence for extended mass 
is seen in the maser rotation curve. The gray stripe represents 
the radial total mass distribution of the best-fit model for the 
CRIRES rotation curves (taking into account errors on Mudm 
and MfL). The orange line represent the total mass distribu- 
tion which is a combination of BH mass, nuclear (4pc size) and 
ex tended mass distribution, with the mass distribution derived 
bv lMueller Sanchez et al. (2006). The nuclear mass distribution 
corresponds to what we measure as Mudm- The green line is 
similar to the orange line but the nuclear mass distribution has 
now a size of 2pc. The dashed green line is the contribution of 
BH mass and nuclear mass distribution, while the dotted line 
represents BH mass and extended mass distribution. The con- 
tinuous black line represents the radial mass distribu tion of the 
nucle ar star cluster at the center of t he Milky Way (|Genzel et all 
120101) combined with the BH mass dGreenhill et alJl2003al) and 
rescaled to match the measurements at ~ 10 pc scale. 



iGreenhill et al.l (l2003al) from modeling of the H2O maser rota- 
tion curve and is located at the inner radius probed by the maser 
rotation curves. The maser disk is almost edge-on. The gray re- 
gion characterized by r < 20 mas is the region probed by the 
H2O maser rotation curve, where the dynamical mass is mostly 
unresolved and belongs to the BH. The outer limit of this re- 
gion is set by the maximum radius probed by the maser rota- 
tion curve. The red filled circle represents the CRIRES spec- 
troastrometry mass measurement obtained assuming a disk with 
i = 61° with respect to the line of sight. This is the inclination 
which is determined from the analysis of the rotation curves at 



larger scales, and is close to the inclination of the galactic disk. 
Since the scales probed by spectroastrometry are intermediate 
between those of the galactic and maser disks, we have also con- 
sidered the possibility that the spectroastrometry is probing ro- 
tation from an edge-on disk (open red circle). Clearly, even if 
the disk were edge-on, the observed kinematics at ~ 2pc-scales 
would still require the presence of extended mass. The black cir- 
cle represents the mass measurement from the CRIRES rotation 
curves, and is located at half the spatial resolution. Finally, the 
brown and orange points represent meas urements of dynamical 
masse s en closed in 8 pc and 9 pc radii bv lMueller Sanchez et al.l 
(120061) and lHicks etatl (120091) . respectively, using SINFONI AO 
assisted observations. It is clear that the measurements from 
SINFONI and CRIRES agree within the errors. The gray stripe 
represents the radial total mass distribution of the best-fit model 
for the CRIRES rotation curves (taking into account errors on 
Mudm and M/L). This is the mass distribution one would have 
derived considering only the K band surface brightness profile. 

Inspection of Fig. QT]reveals the existence of a spatially ex- 
tended mass distribution within the inner lOpc of the nuclear 
region of Circinus. It also shows the capability of spectroastrom- 
etry to improve the spatial resolution by a factor ~ 3, providing 
constraints on the extended mass distribution but, unfortunately, 
it is still not possible to measure the BH mass. Indeed the stellar 
mass distribution recovered from the NICMOS K-band contin- 
uum light profile does not account for the mass distribution in the 
inner lOpc and leads to an overestimated unresolved dynamical 
mass. 

The extended mass within lOpc might be accounted for by 
stars not identified in the NICMOS image due to the presence of 
the hot dust emission or might be accoun ted for by cold gas. The 
AO-as sisted SINFONI observations by iMueller Sanchez et al.l 
(2006) provide the surface brightness profiles of hot dust and 
stellar emission, as well as those of Bry and H2. As shown in 
their figures 3 and 4, the stellar emission is more extended than 
the hot dust and gaseous components. Moreover, the stellar com- 
ponent if offset by ~ ff.'2 and its mass is ~ 1O 6 M making it un- 
likely that Mudm can be accounted for by stars whose emission 
is hidden in NICMOS by the strong hot dust emission. On the 
other hand, they find that Bry and H2 surface brightness distri- 
butions are consisten t with an exponential disk w ith scale radius 
r d =s 4pc. Moreover IMueller Sanchez et ail (120061) find that the 
nuclear H2 luminosity corresponds to a molecular gas mass of 
1.7 x 1O 7 M , fully consistent with the dynamical mass within 
10 pc. Therefore it seems that the extended mass distribution 
within lOpc might be accounted for by a disk of mostly molecu- 
lar gas. In order to verify this hypothesis, we plot with an orange 
line in Fig. fTTIthe mass distri bution obtained by considering the 
BH mass by IGreenhill et al.l d2003al) . the stellar mass distribu- 
tion from the NICMOS K band surface brightness profile and an 
exponential disk with r& = 4 pc and mass 1.7 x 10 7 M H . T his is 
consistent with the analysis of Mu eller Sanchez et al.l (120061) but 
fails to reproduce our spectroastrometric measurements because 
this mass distribution is too extended. A much better agreement 
is found considering an exponential disk with rj = 2 pc and mass 
1.4 x 10 7 M Q which is shown by the green solid line. For com- 
parison, the dashed and dotted green lines represent the relative 
contribution of the gas and stellar masses (each of them added 
to the BH mass). This is perfec tly consistent with the results 
by IMueller Sanchez etaL I (120061) considering the uncertainties 
in converting the emission line surface brightness distributions 
in mass distribution. 

The spectroastrometric measurements probing smaller scales 
thus indicate that the mass distribution is more concentrated 
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than is inferred from the Bry and H2 surface brightness distri- 
butions. The mass which is unresolved even by spectroastrom- 
etry has a large average density of 1.4xlO 5 M pc~ 3 and is thus 
expected to host significant a mounts of star formation. Indeed 
dMueller Sanchez eta found that significant star forma- 

tion is ongoing within the central 8 pc and identified the molecu- 
lar medium with the torus required by the AGN unified model 
and believed to exist in the nucleus of Circinus, a Seyfert 2 
galaxy. 

From the green line it is also possible to infer that the con- 
tribution of gas and stars to the mass distribution equals that of 
the BH mass at 1 pc which is then to be considered as the 
radius of the gravitational sphere of influence of the BH. This 
is perfectly consistent with the independent estimate based on 
the stellar velocity di spersion. For a stellar velocity dispersion 
of cr ^ 75 km/s (e.g., [Graham 2008b) the sphere of influence is 
fBH = GMbh/ct 2 = 1.2pc remarkably similar to what we have 
just found. 

Finally, the black line represents the mass profile which is 
obtained by considering the BH mass and the star cluster in the 
Milky way dGenzel et al.ll2010h wit h a mass reseated to match 
the observed enclosed mass at 1 pc. iGenzel et all d20 1 Oh model 
the Milky way nuclear star cluster with a double power law den- 
sity distribution characterized by a break radius of 0. 25 pc, inner 
power law r~ L3 , outer power law r -18 and density at the break 
radius 1.35 x lO 6 M pc~ 3 . 

The comparison with the Milky Way is suggested by the sim- 
ilarity of morphological types between the two galaxies and in- 
dicates that, overall, the nuclear mass distributions at r < lOpc 
scales are not significantly different. It is intriguing to think that, 
if all the molecular gas in Circinus is converted into stars then 
it will produce a nuclear star cluster similar to that observed in 
the Milky way and also commonly in late type spirals. Since the 
molecular material is likely associated with the torus invoked 
by the AGN unified model, it is also intriguing to speculate that 
there exists an evolutionary track in which dense nuclear molec- 
ular tori are converted into nuclear star clusters. 



8. Conclusions 

We have presented new CRIRES spectroscopic observations of 
the Bry emission line in the nuclear region of the Circinus 
galaxy, obtained with the aim of measuring the BH mass using 
the spectroastrometric technique. The Circinus galaxy is an ideal 
benchmark for the spectroastrometric technique, given its prox- 
imity and the existence of a secure BH measurement obtained 
from observations of the nuclear H2O maser disk. 

The kinematical data have been analyzed using both the clas- 
sical method based on the rotation curves and with a new method 
developed by us and based on spectroastrometry. 

The classical method indicates that the gas disk has inclina- 
tion and position angles consistent with those of the large scale 
galactic disk. The disk rotates in a gravitational potential result- 
ing from an extended stellar mass distribution and a spatially un- 
resolved dynamical mass of (1.7 + 0.2) x 10 7 M Q , concentrated 
within r < 7 pc, corresponding to the seeing-limited resolution 
of the observations. The estimates of inclination, position angle 
and mass are consistent with previous measurements based on 
SINFONI spectroscopy. 

The spectroastrometric method is capable of probing the gas 
rotation at scales which are a factor ~ 3.5 smaller than those set 
by the spatial resolution. The dynamical mass that is spatially 
unresolved at the scales probed by spectroastrometry is a factor 



~ 2 smaller, 1 .9+\\ x 1O 6 M , than that recovered from the ro- 
tation curve analysis, indicating that spectroastrometry has spa- 
tially resolved the nuclear mass distribution down to 2-pc scales. 
This unresolved mass is still a factor ~ 4.5 larger than the BH 
mass measurement obtained with the H2O maser emission indi- 
cating that, even with spectroastrometry, it has not been possible 

to resolve the sphere of influence of the BH. 

Ba sed on the detailed analysis by iMueller Sanchez et al.l 
(2006), this dark mass distribution is likely composed of warm 
molecular gas as revealed by the H2 1-0 S(l) emission at 
2.12yum. This medium has a large density, is clumpy and form- 
ing stars, and has been tentatively identified with the molecular 
torus required by the AGN unified model, which prevents a di- 
rect view of the central BH in Circinus. The mass distribution is 
similar in shape to that of the nuclear star cluster of the Milky 
Way. We tentatively speculate that molecular tori, forming stars 
at a high rate, might be the precursors of the nuclear star clusters 
that are common in late type spirals. 
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